Random cascade model in the limit of infinite integral scale as the exponential of a 
non-stationary 1/f noise. Application to volatility fluctuations in stock markets 



o 



H 
d 



> 
o 



o 



x 



Jean-Frangois MuzjQ 
[/MP ffPj^ CWPS, Umversite de Corse, Quartier Grossetti, 20250 Corte, France and 
CMAP UMR 7641 CNRS, Ecole Poly technique, 91128 Palaiseau, France 

Rachel Bai'ltf] 

SPE UMR 6134 CNRS, Umversite de Corse, Vignola, 20200 Ajaccio, France 

Emmanuel Bacrj{§ 
CMAP UMR 7641 CNRS, Ecole Poly technique, 91128 Palaiseau, France 
(Dated: January 18, 2013) 

In this paper we propose a new model for volatility fluctuations in financial time series. This 
model relies on a non-stationary gaussian process that exhibits aging behavior. It turns out that 
its properties, over any finite time interval, are very close to continuous cascade models. These 
latter models are indeed well known to reproduce faithfully the main stylized facts of financial time 
series. However, it involve a large scale parameter (the so-called "integral scale" where the cascade 
is initiated) that is hard to interpret in finance. Moreover the empirical value of the integral scale is 
in general deeply correlated to the overall length of the sample. This feature is precisely predicted 
by our model that turns out, as illustrated on various examples from daily stock index data, to 
quantitatively reproduce the empirical observations. 
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I. INTRODUCTION 



For several decades, random cascade models have been 
at the heart of a wide number of studies in mathematics 
as well as in applied sciences. They were introduced to 
account for the intermittency phenomenon in fully de- 
veloped turbulence and are involved every time one ob- 
serves a multifractal (or a multiscaling) behavior in the 
variations of statistical properties of some field across 
different scales. Multifractal scaling is generally associ- 
ated with the existence of a random cascade by which 
small scale structures are constructed from the splitting 
of larger ones and multiplication by a random factor. 
One clearly sees that such a scenario necessarily implies 
the existence of a large integral scale T where the cas- 
cade is initiated. As emphasized below (see Appendix 
[A)) , on a general ground, one can show that the moment 
multiscaling behavior of the increments associated with 
any multifractal field cannot hold over an infinite range 
of scales. It necessarily involves a large scale T above 
which scaling properties becomes trivial. In turbulence 
this scale naturally corresponds to the injection scale, 
i.e., the time/space scale where kinetic energy is injected 
into the flow The main question addressed in this 
paper concerns the fields where a multifractal behavior 
is observed without the existence of any obvious integral 
scale. This is notably the case in empirical finance. 
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In quantitative finance, volatility is one of the most im- 
portant risk measures since it corresponds to the condi- 
tional variance associated with price fluctuations at any 
time t @- A well known stylized fact is that volatil- 
ity fluctuations are organized into persistent clusters. A 
huge amount of the econometrics literature is devoted to 
the modeling of this volatility persistence. Among all 
the proposed alternatives, the GARCH models Q and 
their extensions have been thoroughly studied. The ma- 
jor drawback of such models is that, on one hand, their 
aggregation properties are not easy to control and on 
the other hand, they cannot account for the long-range 
nature of volatility correlations d, Q . This last feature 
translates in the fact that GARCH parameters are often 
found to be at the borderline of the stability region. This 
is the so-called IGARCH effect @. 

Under the impetus of early studies of Mandelbrot and 
his collaborators Q , the notions of multifractals and ran- 
dom cascades have been proposed to account for the 
volatility dynamics in many studies of financial time se- 
ries (see e^g. The class of continuous random 
cascades [13[ and in particular the MRW model, provides 
a parsimonious class of random processes that reproduces 
very well most of stylized facts characterizing the price 
return fluctuations 

UBS- Unlike GARCH models, these 
models are continuous time models (so they do not in- 
volve a discrete time step) which aggregation properties 
are easy to handle since they possess some self-similarity 
properties. Within this framework, various empirical es- 
timations reported so far indicate that the value of T 
can vary from few months [1, H| to several years [13 , Il5j 
(see Fig. [H] below). Even if it is well admitted that a 
precise estimation of T can be hardly achieved [3, , 
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one can naturally wonder why one observes such a large 
range in the estimated integral scale values. Beyond the 
problem of the determination of T, a challenging ques- 
tion remains to understand the meaning of the integral 
scale in finance. Unlike turbulence, there is no natural 
large scale that would obviously appear to be associated 
with some "source of volatility" . 

The idea we propose in this paper is that such a 
scale does not exist (or is formally "infinite") and that 
the volatility is a non-stationary process. Let us notice 
that, within standard econometric framework, many au- 
thors already proposed to explain the above mentionned 
IGARCH effect by the non-stationary nature of volatil- 
ity fluctuations: these models include Fractionally inte- 
grated GARC H fTetl . GARCH models with time vary- 
ing parameters [171 fl8l] . stochastic volatility models with 
unit roots [l9j . Our approach is original in the sense 
that we account for the non-stationary nature of volatil- 
ity fluctuations while remaining within the framework of 
multifractal processes. Indeed, we will show that our 
model is such that every single trajectory, for each fi- 
nite time interval, can hardly be distinguished from the 
path of a multifractal process where the integral scale 
is precisely the length of the time interval under consid- 
eration. Our construction is written as the exponential 
of a non-stationary 1/ / noise and is based on an exten- 
sion of continuous random cascades based on infin itely 
scattered random measure as introduced in Refs [13l |20j . 
We show that this process is well defined in the sense of 
distributions and cannot be distinguished from a contin- 
uous cascade model (as the MRW process defined in (Toj ) 
over any finite time interval far from the time origin. We 
check and illustrate our results on some numerical sim- 
ulations. We then consider applications to stock index 
market data that are shown to exhibit some "aging" be- 
havior as precisely predicted by our model. 

The paper is organized as follows: In section [TT| we 
make a brief overview of multifractal models as they have 
been proposed to account for the volatility fluctuations 
in financial time series. The construction of log-infinitely 
divisible continuous random cascades as introduced in 
[l3l [20| is also explained but we mainly focus on the 
log-normal case. In section IIIII we show how one can 
some extend the former cascade models in the formal 
limit T — > +00. The price to pay is that the model is 
no longer stationary. However, this new model has ap- 
pealing properties since, in some sense, it reduces to a 
multifractal model over any bounded time interval with- 
out involving any large scale parameter. Our results are 
illustrated using numerical simulations. In section llVl we 
address the problem of the model estimation using a sin- 
gle realization. We then show that our approach is per- 
tinent to account for the observed volatility correlations 
from intraday to many year time scales. In particular it 
allows one to understand the wide range of estimated in- 
tegral scale values reported in the literature. We use the 
Dow- Jones daily data recorded over several decades to 
provide evidence against the stationarity of the volatility 



process. Concluding remarks and pathes for future re- 
search are provided in section [V] Some technical results 
are reported in Appendices. 



II. MULTIFRACTAL VOLATILITY MODELS: A 
BRIEF OVERVIEW 

A. Multiscaling 

As mentioned in the introductory section, multifrac- 
tal models have provided a family of stochastic processes 
that accounts very well for the main statistical features 
of financial time series [U [22[ . In this section we re- 
call the main results concerning random cascade models 
and set the main notations. We refer the reader to Refs 
[Tol [l3l I20I . [23j for more mathematical details. 

As first proposed by Mandelbrot et al. [7j, multifrac- 
tal processes X(t) with zero mean and stationary in- 
crements, can be constructed through an auxiliary non- 
decreasing multifractal measure M(t) as 

X(t) = B [M{t)\ (1) 

where B(t) is a self-similar process (i.e. such that 
B(Xt) =iaw ^ H B(t)) in general chosen to be a standard 
Brownian motion (H = 1/2). It results that the incre- 
ments of X and M are related by: 

X(t + r)-X(t) = {M{t + T)~M{t)) H X{\) 

law 

= M(r) H X(l) (2) 

Law 

In other words, the variations of M(t) are related to the 
local variance of a Brownian motion. In finance, X(t) 
represents some asset price (or the logarithm of an asset 
price) whose increments are the so-called asset returns. 
In that case, the measure M{t) is usually referred to as 
the "trading time" or the "volatility process" since its 
increments M(t + t) — M(t) > simply correspond to 
the volatility (i.e. the local variance) between times t and 
t + r. Henceforth, most of our considerations will concern 
the "volatility" M(t). All the results can be extended 
to the "price" process X(t) in a straightforward manner 
using Eq. (0). 

In a loose mathematical sense, a non-decreasing 
stochastic process M(t) is called multifractal (or "multi- 
fractal measure") if the moments of its increments (as- 
sumed to be stationary) 5 T M(t) = M(t+T) — M(t) verify 
the multiscaling properties: 

E [5 T M(t) q ] = E [M(r) 9 ] ~ C q r^ , (3) 

where £(q) is a nonlinear concave function of the moment 
order q. Notice that the multifractal nature is prop- 
erly defined by the nonlinearity of ((q) as opposed to 
monofractal situations where is a linear function. In 
order to quantify the multifractality, one often defines the 
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so-called intermittency coefficient A 2 as the curvature of 
£(g) around q = 0: 

A 2 = -C"(0) > . (4) 

The last inequality simply comes from the concavity of 
the spectrum. Indeed, the scaling behavior of Eq. 
([3]) is generally interpreted in the limit of small time 
scales r — > 0. Accordingly, if one computes for exam- 
ple the kurtosis behavior, 

E[M(^] C(4) _ 2C(2) 
E[Af(r) 2 ] 2 

one directly sees that, because J~(t) > 1, one must have 
C(4) < 2C(2). As shown in Appendix El this kind of 
argument can be generalized (thanks to Holder inequal- 
ity) to prove that £(q) is concave. Thanks to Eq. ©, 
one can conclude that the increment probability density 
functions (pdf) of X(t) (the price returns in empirical 
finance) cannot keep a constant shape at different time 
scales r (that would be gaussian in the monofractal situ- 
ation) . It necessarily becomes more and more leptokurtic 
as t — y 0. Both multiscaling and increasing flatness at 
small scales are two well known stylized facts character- 
izing the return time series of many different financial 
markets 0, G3, . 

Let us remark that the previous argument can also 
be used to show that the scaling ([3]) cannot hold for 
arbitrary large scales r. Indeed, since J-(t) > 1, if 
C(4) - 2C(2) < then Eq. (0) can be valid only on a 
bounded range of scales. Therefore there exists an inte- 
gral scale T below which multiscaling holds and beyond 
which one observes monofractal scaling properties (see 
Appendix . 

B. Continuous cascades 

Explicit constructions of multifractal measures can be 
naturally obtained within the framework of random cas- 
cades. The picture of a random cascade comes from the 
physics of turbulence where kinetic energy injected in the 
flow at some large scale is transferred to the finest scales 
by successive steps of eddy fragmentation [l| . The large 
scale where the cascade "starts" corresponds precisely 
to the integral scale introduced previously. Accordingly, 
a discrete random cascade can be constructed as follows: 
one starts with an interval of length T where the measure 
M(dt) is uniformly spread (meaning that the density is 
constant) and splits this interval in two equal parts: On 
each part, the density is multiplied by (positive) i.i.d. 
random factors W. Each of the two sub-intervals is again 
cut in two equal parts and the process is repeated in- 
finitely. Given the discrete and non-stationary nature of 
such constructions and the fact that they are only defined 
in a fixed bounded interval (of size T), more recently, con- 
tinuous cascade constructions have been proposed. These 
models can be viewed as a "densification" of the discrete 



construction [2(| HH, [H| where the multiplication along 
the dyadic tree associated with successive fragmentation 
steps, 

dM = Y[W t = e^ ln ^ , 

i 

is replaced by the exponential of an integral (instead of a 
discrete sum) of a white noise (instead of ln(W)) over a 
cone-like domain in the time-scale plane (instead of the 
tree- node set). More precisely, one defines [l3l[20|: 

dM t , T (t) = Mt, T ([t,t+ dt}) = e" f ' T(t) dt (6) 

with 

Wi,r(*) =Rt+ / dW(u,s) (7) 

where fie,T is a constant such that E [e w<! > r W] = i j 
dW(u, s) is a white noise associated with some infinitely 
divisible law (more precisely an "independently scattered 
random measure" [Tj]) and C^rif) is the cone like do- 
main 1 : 

(u, s) G C e , T (t) {s >£,t- mm(s, T) < u < t} (8) 
This construction is illustrated in Fig. [TJ The final 



s=T 



s=l 




t 1 t 3 

time 



FIG. 1: Construction of a continuous cascade process: uie : T(t) 
is the integral of a white noise over a cone-like domain 
Ci t r(t) in the time-scale plane. The covariance of ui^iti) 
and w<,r(fa) corresponds to the area of the intersection 



Let us remark that the construction we consider here is a "causal 
version" of the original construction proposed in Refs [l3l , |20| 
where a symmetrical cone was used. All the results and compu- 
tations remain unchanged for both constructions. 
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multifractal measure dMx is then obtained as the weak 
limit of dMt T when I — »■ 0, i.e., 

M T (t) = lim / dM iT {t) = lim / e^ T(t) dt (9) 

For the sake of simplicity, we will consider, in this pa- 
per exclusively log-normal random cascades. All our re- 
sults can be easily extended to arbitrary log-infinitely 
divisible laws within the framework introduced in Refs 
[l3|,[2(|. In the log- normal case, dW(t,s) is a 2D Gaus- 
sian (Wiener) white noise of variance X 2 s~ 2 dtds and it 
is easy to see (see Fig. [I]) that the covariance of cje^iti) 
and ui^rih) is simply the area of C^t^i) H C^t^)- Its 
expression reads: 

f A 2 ln(^) if£ < t < T 
Cov (u t>T {t), ut,T(t + t)) = \ A 2 (ln(f ) if r < £ 

[ if r > T 

(10) 

In that respect the mean value of lu^t has to be chosen 
as: 

A 2 / T \ 

W,T = -y f l + ln(-)J . (11) 

Notice that in the log-normal case, the intermittency co- 
efficient A 2 and the integral scale T are the only parame- 
ters that govern the multifractal statistics. The previous 
equation mainly says that the logarithm of a random log- 
normal multifractal measure is a Gaussian process which 
covariance decreases as a logarithmic function, log(T/r). 
This features has been shown to directly reflect the tree 
structure of discrete random cascades (see Refs [9l [26|). 



C. Stochastic self-similarity 



One can establish another self-similarity property [27[ 
when one also rescales the integral scale. In that case, 
one has trivially from Eq. (fTU|) . V r > 0: 

U>rt,rT(rt) = Ui >T (t) (15) 
law 

M rT (rt) = rM(t) (16) 

law 

which means that a trivial scaling is obtained when the 
integral time T is rescaled with the time. 

In the field of empirical finance, random cascades have 
allowed one to understand that the observed multiscaling 
properties of return moments and the long-range corre- 
lated nature of the volatility are the two faces of the same 
coin. The (log-normal) multifractal random walk model 
has proven to be a simple, parcimonious model that re- 
produces most of observed statistical properties of asset 

returns H OH Hum. As far 

as statistical estimation is- 
sues are concerned, as shown in Rcf. [291 ] . intermittency 
exponent estimations based on Eq. (|10[) are far more re- 
liable than those based on moment multiscaling ([3]) (see 
also [HI, H3| for additionnal results on the intermittency 
exponent estimation using GMM methods). Empiri- 
cal evidence for the logarithmic nature of log- volatility 
correlations have been provided for different asset price 
time series over different markets [H, [HI EH [H|. 
All these results confirm the multifractal nature of as- 
set return fluctuations with an intermittency coefficient 
A 2 G [0.01,0.03]. However the reported values of the in- 
tegral scale T vary in a wide range of scales, between few 
months and several years. The main question we want 
to address in this paper concerns that point: what is the 
value of the integral scale in financial time series ? 

III. THE LIMIT OF INFINITE INTEGRAL 
SCALE: A NON-STATIONARY MODEL FOR 
LOG- VOLATILITY 



All the (multi-)scaling properties of M(t) (and subse- 
quently of X(t)) can be shown to result from the loga- 
rithmic nature of this covariance. Indeed, since u^t(*) is 
a Gaussian process, one can directly infer from Eqs. (fT0|) 
and HT]) that, V r < 1, V t < T, 

IJ->r£,T(rt) = U)£ t T(t) + (12) 
law 

where Q r is a normal random variable of variance 
—A 2 ln(r) and mean 4r ln(r). From Eq . (fT2l . the stochas- 
tic self-similarity property results [13|, [2(| : 

M T {rt) = re nr M{t) (13) 

law 

which directly proves the multiscaling (Eq. (|3])) of the 
moments of M(t) (and thus of X(t)) with a parabolic 
function: 

i nE r e 9n.i x 2 AV 

C(?) = 9+^ r ^ = ff (l + T )-^-. (14) 



A. Definition of the model 

The broad range of observed values of the integral 
scale in empirical studies leads us to ask the question 
of the interpretation of the integral scale value in finan- 
cial markets. Unlike turbulence, there is no obvious large 
scale that could be singularized and associated with some 
"source" of randomness. Even if the heterogeneity of 
agents and the wide range of time horizons used by mar- 
ket participants is a well recognized fact, this can hardly 
be invoked to define a single scale that could be as large 
as several years. 

A way to answer the previous remarks could be to con- 
sider the model introduced in (3l| where the authors re- 
placed the log-correlated uie,T(t) by a long-range (e.g. a 
fGn) correlated stationary Gaussian process. However 
the continuous time limit of such a process is trivial (i.e., 
it necessarily involves a small scale cut-off) and its scal- 
ing properties are not exact and hard to handle. Another 
solution is to define a random cascade process in the limit 
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T — > oo. However, the definition of such a limit is not 
obvious, since, as emphasized in the previous section and 
shown in Appendix El one cannot define any multiscaling 
behavior without involving a finite integral scale. As one 
can check in Eq. (fTOj) . by letting T —> oo, one obtains 
an infinite value of the variance (and the mean) of toe,T- 
In Ref. [lH, the authors have considered the possibil- 
ity of an infinite integral scale and provided an explicit 
prediction formula of uji.t^oo (that we denote as uJe.oc). 
However this process is not defined in a classical sense 
but only in some quotient space, namely a space of pro- 
cesses defined up to constant time functions. It has been 
shown that 

lim / (j>{u) oji x(t — u)du (17) 

T^oo J 

is meaningful for a class of smooth functions (f> provided 
it is of zero mean. We already know that the singularity 
of the covariance function at r = when £ — » (Eq. 
(fTU)) ) means that the limit of w^t (or exp(w^T)) has to 
be considered as a noise process and is well defined only 
when interpreted in the weak (distribution) sense. When 
T — > +oo, Duchon et al. [28| show that u^_>o,oo can be 
still interpreted in a weak sense but only for test functions 
satisfying J (j>(t)dt = 0. This process and notably its 
exponential e" f '°°, is however hard to interpret and of 
unclear practical interest in quantitative finance. 

In order to handle the low-frequency problem related 
to T — > +oo, we propose in this paper an alternative solu- 
tion that consists in considering a non- stationary process 
where, at time t, the integral scale is precisely T — t. We 
define a process u)i{t) as for standard cascade, from the 
integration over a cone-like domain in a time-scale plane, 
where, at time t, the parameter T in Eq. ([8]), is replaced 
by t; 

(u,s)€C<t(t) O {s > £,max(0,t - s) < u < t} (18) 
if t > I 

C e (t) = otherwise . (19) 
The process u>e(t) is then defined by: 

u e (t) =in(t)+ [ dW{u,s) (20) 

J(u,s)eC t (t) 

where fig(t) is a deterministic mean value defined be- 
low and dW(u, s) a Gaussian white noise of variance 
\ 2 s~ 2 duds. This construction is illustrated in Fig. [5] 

In Fig. [3] we have plotted a sample of uii(t) generated 
at rate At = £ = 1 over 500 points. As one can see, the 
non-stationary nature of w^(t) is not obvious (see below). 

We can compute the covariance of u>i(t) that corre- 
sponds to the domain Q intersection areas (see Fig. [2]). 
For t\ < t% = t\ + t, its expression reads: 

( A 2 ln(^-) if t > I 

Cav(Mti),Mt2)) = { x2 Nt) + 1 - f) if t < ^ 
[ \it<l 

(21) 




time 



FIG. 2: Construction of the non-stationary UJ((i) process as 
the integral of a white noise over a cone- like domain Ct(t) in 
the time-scale plane. The covariance of LJi(ti) and uefa) is 
simply the area of the intersection C'e(ti) PI C(,(t-z) 




200 400 

t 



FIG. 3: An example of a path of UJi(t) where the numerical 
construction has been performed by sampling both space and 
scale in Eq. ((20)) . 

This equation implies notably that 

Var[^(t)] = A 2 (l + ln(^)J . (22) 

One clearly sees that (Jg(t) is a non-stationary gaus- 
sian process but its covariance has striking similarities 
with the stationary situation (Eq. ([TU| ) where the in- 
tegral scale has been replaced by the current time t (or 
max(ii,t2) in the covariance expression). 

The non-stationary behavior of the covariance is illus- 
trated in Fig. |31 where we have plotted the estimated 
as well as analytical Cov (u>£ (£i), u>£ (£2)) as a function of 
ln(r) = ln(i2 — ii) for different values of ti- Remark that 
this kind of non-stationarity is reminiscent of an aging 
behavior as observed in off-equilibrium relaxing systems 
[32l [HI where the "age" of the process t^ controls the 
characteristic correlation length. The logarithmic behav- 
ior of the variance is illustrated in Fig. [5] 

Let us show that one can choose a function /J,((t) in 
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> 
o 
u 




provided, V t, 



E 



a condition equivalent, in the log-normal case, to 



1 



A 2 



-Var[w/(t)] =-— (1 + Hj) 



,t 



(24) 



In Appendix [B] we provide an alternative direct proof of 
mean square convergence. 

Notice that this equation also guarantees that 



E[M(t)} = Yax[X(t)] = cr 2 t 



(25) 



(recall that X(t) = B[M{t)\ with B{t) a standard Brow- 
nian motion). 



FIG. 4: Covariance function Cov (uje(ti), oj( (£2)) as a function 
of ln(r), with t = |t a - ti\ and t 2 = 10,40,150,500 (from 
bottom to top curves). The bold lines correspond to numer- 
ical estimates using 500 samples of u)e(t) while the thin lines 
correspond to the analytical expressions (Eq. (|2f [) )■ We have 
chosen I — I and A 2 = I . 





ln(t) 



FIG. 5: Variance Vaxwe(t) as a function of t (a) and \n(t) (b). 
We have superimposed to the expected analytical expressions 
the estimated variance using 500 Monte-Carlo samples 



of cJe(t) with I = I and A 



Eq. PD|) such that one can define the limit £ — > of e u 
in the weak sense, i.e.,: 



M(t) = lim / 
Jo 



(23) 



In fact, as for continuous stationary cascades, one can 
use a general argument on positive martingales (as e.g., 
in Ref. [13]) if, for all time interval /, fje^dt is a 
martingale as a function of I. This is precisely the case 



B. Scaling and self-similarity properties 

Let us remark that increments of we(t), Sh^e{t) = 
u>e(t + h) — uii(t) (h > £) have a time dependent vari- 
ance so they are not stationary. However, for r > h, 
their covariance depends only on the lag r. After a little 
algebra, their expression reads: 



(26) 



Cov{S h ui(t),S h ui(t + T)) = A In 1 



This covariance function corresponds to a power- 
spectrum such that Ps hU (f) ~ l/l when / <C h~ l . Since 
Put(f) ~ f~ 2p 6 h u e (f)> it; res ults that lim^ uj e (t) can be 
associated with a 1// power-spectrum. Let us mention 
that, in Ref. [33j, the author has already raised the pos- 
sibility of an "aging" non-stationary model in order to 
handle the low-frequency problem of 1/f noise. In that 
respect, ui^oit) can be interpreted as the limit H — > 
of a fractional Brownian motion (fBm) Bjj{t) of Hurst 
parameter H [34j . 

This interpretation of uji(t) can also be suggested from 
its self-similarity properties. Indeed, from the covariance 
expression (|2"Tj) , one can establish the following invariance 
relationship for Lde(t): 



(27) 



This equality extends to H = the standard fBm 
self-similarity Bn(rt) =i aw r H BH{t) 2 . It is notewor- 
thy that u!e^o(t) has the drawbacks of both fractional 
Gaussian noises and fractional Brownian motions since 
it exists only in the sense of distributions and it is a 
non-stationary process. From the definition (|23p and 



2 Remark that, since fBm's are defined from both their self- 
similarity properties and the stationarity of their increments ; 34] , 
in full rigor ui£^o(t) cannot be identified to Bu(t) with H = 0. 
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thanks to previous equality, one deduces the simple self- 
similarity property of the volatility process M(t): 



M{rt) = rM{t) . 

law 



(28) 



Let us remark that relation ([27]) is different from Eq. 
(]T2"]) but can be understood as reminiscent of Eq. ([15)) 
where one allows the integral scale to become infinite 
(i.e., T —> oo). 

When one compares the self-similarity of M and Mt 
(Eqs. (|2"8|) and (]T3"]) L one sees that in the former case 
there is no stochastic factor e° r and the scaling of the 
moments of M (and therefore of the return process X(t)) 
becomes trivial: 



E [M{r) q ] = C q r q E [\X{r)\ q ] = K qT i 



(29) 



In the sense of Eq. ©, it thus appears that M(t) (or 
X(t)) is not a multifractal process. However, one must 
carefully interpret the previous equation since M{t) (and 
then X(t)) has no stationary increments. It results 
that there is no reason that the moments E [M(r) q ] and 
E [[M(t + t) - M{t)] q ] behave in the same way. Let us 
make the explicit computation for q = 2. In that case, 



E [M(r) 2 



lim 

£->0 



lim 

t-yO 



dudv 



max(w, v) 



n jo \ \u-v\ 
1 '" / max(u,v) 



C 2 T 



JO 
2 



\u — v\ 



dudv 
x 2 

dudv 



whereas, 



E[[M(t + T )~M(t)} 2 ] = 
lim^o Jt +T It +T E [e^+^] dudv 

-irir {^j^ffdudv 

= * 2 UU{ 1± l^r 1 fdudv. 

If one supposes that j < 1, then in the last integral 
the term max(u, can be neglected and, using the 

change of variables v! — ut/r and v' = vt/r, one gets: 



E[[M(* + t) -M(t)} 2 ] ~ C 2 (t)r 



2-y 



where the constant C^it) ~ t . The previous equa- 
tion shows that in the limit r<Ci, the mean square of 
the increments of M(t) behaves as the increment of the 
multifractal measure Mx(t) (with the scaling exponent 



C(2) = 2 — A 2 ) where t plays precisely the role of the 
integral scale T. This behavior can be directly estab- 
lished from the expression of the covariance, Eq. (|21|) : 
Indeed, let us consider two times t\,t2 in some interval 
[t , t + At]. If At < t , then to the first order in t /At, 
we have Cov (u>i(ti),U)i(t2)) — A 2 ln(io/|ti — ^|), i-e. the 
same covariance as the process u>e,T used to build an exact 
multifractal random measure with T = to- This means 
that the non-stationary process M(t) defined in Eq. (|23p . 
cannot be distinguished from a (stationary) multifractal 
random measure M to (t) of integral scale T = t Q over any 
interval [to, to + At] (to the first order in to/ At). 



IV. APPLICATION TO FINANCIAL DATA 

As recalled in the introduction, various authors have 
suggested that most of stylized facts characterizing the 
volatility associated with asset prices in financial mar- 
kets can be accounted by multifractal measures. Let us 
illustrate how the model M(t) introduced in this paper, 
allows one to explain the large discrepancies of the re- 
ported integral scale values as a consequence of the non- 
stationary nature of log-volatility. Since the model is 
non-stationary and since in practice there is no possibil- 
ity to have an ensemble of many independent samples, 
one has first to discuss which kind of estimation one can 
perform on a single realization of the volatility. 



A. Pathwise properties and estimation issues 

Let us suppose that one studies a multifractal random 
measure Mx(t) (i.e. a classical random cascade with fi- 
nite integral scale T) over an interval [to , to + At] (or, 
since Mt has stationary increments, over [0, At]) with 
At < T. Then from the self-similarity relations (TT21 and 
([15]). one as, for all r < At/T < 1, 



M T (t) = r^Mrrirt) = r^e^ M rT (t) 

law law 



(30) 



Since the random variable uj r is fixed on a single realiza- 
tion, this equality clearly means that one cannot distin- 
guish over any interval [to , to + At] two multifractal mea- 
sures M Tl (t) and M Ta (t) with T x ^ T 2 and Ti,T 2 > At. 
Estimating the integral scale on a single realization of 
Mx{t) over an interval of length At < T is thus im- 
possible. The question is to which value an empirical 
estimation leads to ? 

Empirically, as advocated e.g., in Ref. [14]], the cor- 
relation properties of LOi^{t) can be estimated using a 
proxy (called the "magnitude process") of uJh,T{t) esti- 
mated from the logarithm of the increments of Afr(t): 
uih.T — hiShMxit). If uJh,T is sampled at rate h over a 
time period of length At, the estimator of its covariance 
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Ca*(t) at lag r = nh, reads: 

CA t (r) = (N-n)- 1 (u) h , T [ih]-pL){u h>T [{i+n)h]-fi) 

i=0 

(31) 

where iV = is the sample size and fi is the empirical 

mean: ft = TV -1 J2k=o LjJ h,T{kh). In Appendix ICl fsee 
also Ref. [l4[ for a more technical approach) it is shown 
that: 

This equation means that, over a sample of size At, the 
estimated auto-covariance of the magnitude associated 
with a multifractal process of integral scale T > At is 
the auto-covariance of a multifractal process of integral 
scale e _3/2 At. 

If we now go back to the non-stationary process M(t), 
since we have shown that, over every interval [to, to + At], 
M(t) =iaw Mt Q {t), we can conclude that, as soon as 
to > At, the estimated auto-covariance of LUh(t) = 
\n[M(t + h) - M(t)] will be provided by Eq. (J^J). In 
other words, for observations far from the time origin, 
the estimated integral scale is always (up to a constant 
factor) the overall sample size. This is illustrated in Fig. 
IHIc) where we have reported the estimation of the mag- 
nitude auto-covariance for various sample lengths At. 
More precisely, we have generated a single large sample 
of the process M(t) from which the magnitude time se- 
ries u)h{t) has been computed. This series (of overall size 
L = 2.10 4 ) is displayed in Fig. [UJa). For each subinterval 
size At = 16, 32, . . . , 512, the sample is divided in L/ At 
sub-samples of length At. The reported estimator C&tij) 
is the average of the obtained empirical covariances over 
all of the Lj At intervals. One can check in Fig. [SJc) 
that the theoretical predictions (I3"2l (solid lines) are, for 
all At, in good agreement with the observations (•) and 
one clearly observes an apparent integral scale that grows 
with At (as e" 3/2 At). 

B. Application to daily stock data 

Let us now apply the previous analysis to real data. 
We report below the empirical results we obtained on 
three stock indices (namely the Dow- Jones, the CAC40 
and the FTSE100 indices) over sufficient long time pe- 
riods. In each case, h = 1 day and u>h(k) at day k is 
estimated as LOh(k) = ln(cr(fc)), where cr(fc) is the relative 
range computed from highest and lowest stock values ob- 
served during the day k. The considered time periods 
are 1929-2011 for the Dow- Jones series (around 21.000 
trading days), 1990 to 2011 for the CAC40 (around 5500 
trading days), 1984 to 2011 (around 7000 trading days) 
for the FTSE100. 

In Fig. |UJb) is plotted the time series corresponding 
to the daily log- volatility Ldh(k) of the Dow- Jones index. 
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FIG. 6: (a) Sample path of u>h{t) of length 2.10 4 where the nu- 
merical construction has been performed by sampling in both 
space and scale the cone-like sets Ct(t) for small I in order 
to define M(t). (b) Magnitude estimated as ln(cr(t)) where 
er(t) is the daily range (difference between highest and lowest 
daily return values) associated with the Dow- Jones time series 
from 1929 to 2011. In (c) and (d) are reported the magnitude 
auto-covariance estimation in semi-logarithmic scale for dif- 
ferent sample sizes At = 16, 32, 64, 128, 256 and 512 (from 
bottom to top), (c) Estimation from the model sample path 
in (a), (d) Estimation from the Dow- Jones daily data in (b). 
The solid lines represent the theoretical predictions from Eq. 
[[32} with A 2 = 0.01. 



Very much like the model (Fig. [Ha)), one can observe 
excursions away from the mean value lasting for several 
years. For each of the 3 volatility series, we reproduced 
the same covariance estimation experiment we conducted 
for the model (Fig. [BJc)). In Fig. [^d) are reported the 
results obtained for the Dow- Jones index while in Fig 
[7] are reported the results obtained for the CAC40 and 
the FTSE100 time series. Tough these latter series have 
a smaller size and lead to more noisy results, it clearly 
appears in all cases that the empirical auto-covariance 
functions are fairly well fitted by a multifractal loga- 
rithmic shape A 2 ln(T/r) (i.e. they are linear as func- 
tions of ln(r)) with a constant intermittency coefficient 
A 2 ~ 0.01. However the apparent integral scale T (the 
intercept of each curve) appears to strongly depend on 
At. In Fig. [HJd), we see that the model predictions (solid 
lines) as described by Eq. ([32| . fit fairly well the data. 
This is confirmed in Fig. where we have plotted (in 
log-log scale) the estimated integral scale as a function 
of the sample size At. One can see that the analyti- 
cal prediction T(At) = e~ 3/2 At (solid line) is m very 
good agreement with the empirical data. These results 
allow us to understand the origin of the wide range of 
integral scale values (from few months to several years) 
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reported in the literature so far. This is illustrated in 
Fig. IH1 where we have reported the estimated values of 
the integral scale T gathered from the recent literature 
[8l4l(il [Til |2ll l35j . Even if these studies concern various 
data sets at different time resolutions (intradays, daily,..), 
different time periods and correspond to different asset 
classes (FX rates, stocks,..), we see that the reported val- 
ues of T are spread closely around the theoretical curve 
(solid line in Fig. [9]). 




2 4 2 

ln(-r) ln(-r) 



FIG. 7: Magnitude covariance estimation C At (t~ ) as a func- 
tion of ln(r) for respectively (a) CAC40 (5500 data points) 
and (b) FTSE 100 (7000 data points) stock index series. Each 
curve corresponds to a different sample size At used for the 
estimation (At = 16, 32, 64, 128, 256) and for each At, (5a* (t) 
has been obtained as the mean value over all available peri- 
ods of size At. One clearly sees that the behavior is the same 
than for the Dow- Jones index in Fig. Od): the integral scale 
(intercept) is growing as a function of At. The noise ampli- 
tude is greater because the overall sample sizes are smaller 
than for the Dow- Jones series. 




4 ln(At) 6 



FIG. 8: Estimated integral scale T as a function of the sample 
size At (in log-log coordinates) for the Dow- Jones daily time 
series from 1929 to 2011. The solid line represents the value 
e -3 / 2 Ai one expects theoretically. The reported error bars 
correspond to standard deviation of the empirical mean values 
estimated from the observed dispersion over all sub-intervals. 



V. CONCLUSION AND PROSPECTS 

To conclude we have introduced a new model of 
stochastic measure as the exponential of a non-stationary 




6 8 10 

ln(At) (day) 



FIG. 9: Estimated integral scale T as a function of the sample 
length At (in log-log coordinates) gathered from several recent 
studies in the literature: (•) from [sj, (o) from l2lll . (A) from 
0, (a) from [H, (□) from [l(| and (■) from [lj. The solid 
line represents the model fit ln(T) = ln(At) — 3/2. 

gaussian 1// noise. We have shown that, over any fi- 
nite time interval, provided the considered time t is large 
enough, this model can be hardly distinguished from a 
multifractal random cascade with an integral scale that is 
equal to the sample length. Our approach can be very ap- 
pealing to model all phenomena where multiscaling prop- 
erties are observed without the existence of any natural 
large "correlation" (or "injection") scale in space or time. 
For example, in finance, the agreement of the model pre- 
dictions with the observed behavior of log- volatility cor- 
relation in various stock indices is striking. These find- 
ings suggest a peculiar (aging) non-stationary nature of 
volatility fluctuations. The question of the meaning of 
the time origin, the possibility to estimate this time from 
empirical data will have to be considered in future works. 
On a more general ground, the explanation of such non- 
stationarity is an important question that will have to be 
addressed from the market dynamical properties at mi- 
crostructure level but also within the framework of agent 
based approaches including behavioral finance or theory 
of self-referencing dynamics of market prices. On a math- 
ematical ground, it will be interesting to study this model 
and its possible variants in relationship with fractional 
Brownian motion, since it offers the possibility to give a 
meaning to the limit H — > 0. Finally our approach can 
shed a new light in the field of 1/ f noise modeling. 



Appendix A: Proof of the concavity of f (q) and the 
existence of an integral scale 

Let us prove that if Eq. ([3]) holds in the limit of small 
time scales r, then i) £(q) is a concave function of q and 
ii) it necessarily involves a bounded scale T below which 
it can no longer be valid. We start by assuming that the 
following scaling holds in some range of scales: 

E[|Mr(f)|«]~c g r«d 
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Let F(q, r) = In (E [\S T X(t)\^]). Then by Holder inequal- 
ity, F(q, t) is, for each r, a convex function of q. If one 
assumes it is regular enough so that its second derivative 
exists, one thus has, V r > 0: 



F"(q,T) > 



(Al) 



If the previous scaling law holds, this can be written as: 



d 2 lnC g 
dq 2 



+ C(q)Hr)>0 



(A2) 



One sees that if the scaling holds in the limit r — > 0, this 
inequality can be true only if £"((/) < 0, i.e., must be 
a concave function of q. If the scaling is also valid at scale 
t = 1 (up to a redefinition of r we can always assume 
it is the case), C q is the order q moment of the random 

variable 5iX(t) and c q — d ]?2 q > 0. This means that: 



ln(r) < 



(A3) 



In other words, if Q(q) is strictly concave (multifractal 
case), the scaling can only hold in a limited range of 
scales and there exists an integral scale 



T = inf | e<"w 

9 



above which it is not valid. 



Appendix B: Proof of the mean-square convergence 

of Mi(t) 

Let us provide a direct proof of the mean square weak 
convergence of Mg(t) (or Mi(I) = J T e U( ^du for a any 
given time interval /) as defined in (|23[) when £ — > 0. For 
that purpose let us show that 

\xmE[{M t {t) - M e [t)f) = 0. (Bl) 

Without loss of generality, we assume in the sequel that 
£' > £. Since, 



E[(M e (t)-M e ,(t)) 2 ] =E 



(/o ( eW ' (u) ~ e^' (u) ) du 



where we denoted Cej'(u,v) = Cov (ujg(u), ujfj (v)) and 
used the obvious property Cg t e> (u, v) — Ci>ji (u, v) if £' > 
£. Let us split the integral in 3 domains: 



// + // 

J Je<\u-v\<e' J J\u- 



>o Jo J J\u-v\<e J Je<\u-v\<e' J J\u-v\>e' 
It is clear that in the last interval, Cg,^{u,v) = 

Ci'^'(u,v) = X 2 hi( m | i *^"j"- > ). The corresponding inte- 
gral in Eq. (|B2[) is thus zero. In interval I < \u — v\ < £', 
thanks to expression ([2~Tj) . one has 

( e C ltl (u,v) _ e c e ,^,(u.v)\ dudv = (^i-a 2 ) 

n<\u-v\<t' ^ ' 

while in the last interval, 



z C e , e (u,v) _ e C,/ i( /(«,«) _ Q^ 1 -* 2 *) 



l\u-v\<£ 

proving the mean square convergence (|B1[) . 

Appendix C: Magnitude covariance estimation 



Let us establish Eq. (|3"2"1) . Let us denote C(n) = 
CAt(T = nh) and C(n) — A 2 \n(T/nh) the theoretical co- 
variance as given by Eq. (|10j) at lag r = nh. By taking 
the expectation of expression pip after expanding the 
expression of /t, one finds: 



E 



C{n) = C{n) + K(0) - 2K{n) 



(CI) 



where 



N-n—l N-l 



K(n) 



N(N -n) 

y ' i=Q j=0 



E E c (N-ii) 



If h is small enough, one can replace the double sum by 
its integral approximation: 



K(n = r/h) 



X 2 



At 2 (l - Jo 



At /-At-7 



du dv In 



T 



\u — v\ 



= JIjI du dv E [ e m(u)+w e (v) + e w f /(«)+w < /(«)_ 2e wK«)+"/'(«) 

and since (ujg, ujpj ) is a vector of correlated Gaussian pro- 
cesses, thanks to Eq. (|MJ). E [(M e (t) - M t {t)) 2 ] reduces 
to: 



Evaluating this integral leads, to the first order in r/At, 
to the expression: 



/Tp 3 / 2 
K(n = r/h) = A 2 ( In [ 



2At J 



o 



du dv 



C t ,i(u,v) _ C e , e ,(u,v) 



(B2) 



V At 

Inserting this expression in Eq. (jClj) . one gets Eq. (f32j) . 
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